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o\ : 

^ An analytical theory is presented for the damping of low-frequency adsorbate vibrations via 

resonant coupling to the substrate phonons. The system is treated classically, with the substrate 
^ modeled as a semi-infinite elastic continuum and the adsorbate overlayer modeled as an array of 

^ , point masses connected to the surface by harmonic springs. The theory provides a simple expression 

■ for the relaxation rate in terms of fundamental parameters of the system: 7 — miJjQ / AcpCT , where 

m is the adsorbate mass, ljq is the measured frequency, Ac is the overlayer unit-cell area, and p 
and ct are the substrate mass density and transverse speed of sound, respectively. This expression 
is strongly coverage dependent, and predicts relaxation rates in excellent quantitative agreement 
with available experiments. For a half-monolayer of carbon monoxide on the copper (100) surface, 
the predicted damping rate of in-plane frustrated translations is 0.50 x 10^^ s~^, as compared to 
0^ ' the experimental value of (0.43 ± 0.07) x lO^'^ s~^. Furthermore it is shown that, for all coverages 

presently accessible to experiment, adsorbate motions exhibit collective effects which cannot be 
treated as stemming from isolated oscillators. 
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^2 ' Low-frequency adsorbate vibrations have important consequences in many fundamental chemical processes at sur- 
faces, including reactivity, surface diffusion, and desorption. Understanding the dynamical nature of the low-frequency 
modes can provide important insights into the physical behavior of these processes. When the frequency of an adsor- 
bate mode is so low that it lies within the phonon band of the substrate, it can resonantly mix with the phonons to 
form new normal modes with composite adsorbate- and bulk-like character. Thus the adsorbate mode can be viewed 
^ as developing a finite lifetime by decaying into propagating bulk elastic waves. This resonant coupling is an important 
mechanism in understanding energy relaxation processes at surfaces. It has been shown, for example, to dominate 
^ ' the relaxation dyaamics of the in-plane frustrated translational motion of a carbon monoxide adlayer on the (100) 
' surface of copper.El 

\^ In this article, we present an analytical theory of the resonant coupling of an adsorbate overlayer to its substrate. 

This theory results in simple mathematical expressions for the vibrational damping rate and the renormalized vi- 
brational frequency of the adlayer, in terms of fundamental properties of the system, such as the mass density and 
• speed of sound of the substrate, the mass of the adsorbate molecules, and the adsorbate vibrational frequency. When 
applied to systems for which there are experimental data, these expressions give quantitatively accurate results, thus 
verifying their predictive power. 

In this theory, we model the substrate as a flat, isotropic, semi-infinite, classical elastic continuum. The elastic- 
"j^ ' continuum approximation is valid when considering adsorbate modes lying in the frequency range of long-wavelength 
I substrate acoustic phonons, where the dispersion can be accurately modeled with a Debye spectrum. Furthermore, 
I the choice of a flat and isotropic substrate is particularly appropriate for low-index surfaces of close-packed materials. 
''O Thus, this theory is especially applicable to the in-plane frustrated-translational modes of molecules on low-index metal 
Ch , surfaces.EJ These modes typically have frequencies in the acoustic range of the substrate. The molecular overlayer is 
O • modeled as a two-dimensional periodic lattice (or net) of identical, classical point masses, each connected to the 
surface of the substrate by a classical harmonic spring. p. 

The theory presented here builds on the earlier work of Persson and Ryberg (PR) ,0 who examined resonant vibra- 
tional decay for the case of a single, isolated adsorbate coupled to an elastic substrate. They deduced the following 
expression for the damping rate: 
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7=--p cooC, (1) 

where loq is the resonance frequency, m is the adsorbate mass, p is the substrate mass density, ct is the transverse 
speed of sound in the substrate, and ^ is a dimensionless constant of order unity which depends only on the ratio of 
the transverse and longitudinal speeds of sound. We show below, however, that a dense layer of adsorbates behaves 
fundamentally differently. Collective motion of the adlayer strongly enhances the coupling efflciency compared to 
the isolated case, whereas interference effects restrict the phase space of substrate modes available for coupling. One 
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might imagine that both hmits could be achieved experimentally by varying the surface coverage. However, we also 
show below that the dense-adlayer limit holds for an enormous range of coverages. The isolated-adsorbate limit does 
not set in until coverages well below that which can be observed by current experimental techniques. 

We begin by developing the theory in its most general form, which includes both the dense-adlayer and the isolated- 
adsorbate extremes as limiting cases. Afterwards, we indicate how to deduce the PR result and then focus in detail 
on the dense-adlayer limit, as this is the experimentally realizable case. In a forthcoming paper, we will address the 
general result in greater detail, examining the cross-over behavior between the two extremes. 

Consider a semi-infinite, classical elastic medium described by the displacement field u(x, t) and occupying the half- 
space z > 0. Its surface is the xy-plane. (In this paper three-dimensional vectors are denoted in boldface, whereas 
two-dimensional vectors on the surface are designated with arrows.) The surface is covered by a two-dimensional 
periodic array of point masses m attached at lattice sites {Ri} by classical harmonic springs of bare frequency ujq 
(i.e., spring constant majp). Let Si{t) denote the displacement coordinate of the ith adsorbate from its equihbrium 
position. In general, the polarization of Si can have both in-plane and normal components. However, for simplicity, we 
consider Si to be fully a;-polarized. It can be shown, in fact, that the mathematical formulas describing the dynamical 
response of in-plane and normal motion are identical up to numerical constants of order unity. Furthermore, x- and 
y-polarized motion are degenerate by symmetry. 

The equation of motion for the ith adsorbate is given (in the frequency domain) by: 

-muj'^Si{LL>) + mujl s^{uj) ~Ux{Ri,uj) = /^(w) , (2) 

where /i(a;) is the a;-component of an external driving force on the ith adsorbate. If the substrate were rigid (i.e., 
u(x,u)) = 0), the adsorbates would clearly undergo simple harmonic motion at frequency loq. However an elastic 
substrate couples to the adsorbate motion, thereby damping the vibration and renormalizing its frequency. The 
influence of the coupling on the adsorbate dynamics can be described mathematically by the complex poles of the 
adsorbate response function, x{RiT^)i defined by 

Si{i^) = x{Ri,^)ft{^) ■ (3) 

To determine x(^i,w), we first integrate out the substrate degrees of freedom, Urc{x,u}). _ 
The equation of motion for an elastic substrate subject to an external driving force, a{x,uj), at the surface isu 

- Lj'^u{x,uj) - c|V (V ■ u{x,uj)) + CyV X (V X u{x,uj)) = 6{z)a{x,uj)/p , (4) 

where cl and Ct are the longitudinal and transverse speeds of sound, respectively, specifying the elastic properties of 
the material. The surface driving force is provided by the oscillating adsorbates, and thus is given by 

ax{x,uj) = y^TOtjQ Si{uj) - u^,{Ri,uj) 5{x - Ri) (5) 

i 

ay^^{x,uj) = . 
The solution to Eqs. (|) and (||) is given formally by 

Ua(x,z]Lo) = I Dapix, z; x' , z';ll>) af3{x' , z']uj) (Px'dz' , (6) 

where a and (3 run over Cartesian directions, and Dap{x, z]x' ^ z' ]uj) is the Green's function of the elastic substrate, 
defined as the a-component of the solution to Eq. (jj) with source term 5'^{x — x')6{z — z')/p in the /3-direction and 
stress-free boundary conditions. Several simplifications can be made for the present problem. First, the driving force, 
a, exists only at the surface and is x-polarized. Furthermore, according to Eq. (|^), only the x-component of u{x,ll!) 
at the surface is of interest. Finally, translational symmetry parallel to the surface means that I?a/3 can only depend 
on the difference x ~ x' . Thus, Eq. (0) becomes 



Ux{x,uj) = j Dx.j;{x - x\uj) ax{x' ,uj)d x' , (7) 

where the z and z' arguments are zero and have been suppressed for convenience. In reciprocal space, 

u^{k,uj) = D^.j,{k,uj)ax{k,uj) . (8) 
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Because of the two-dimensional periodicity of the adsorbate lattice, the stationary solutions to both the substrate 
and adsorbate equations of motion must obey Bloch's Theorem. Namely, 



Si{u) = sq{uj) e"^'^' , 

where g is a wavevector in the first Brillouin zone of the adsorbate lattice. Combining Eqs. (^) and 
transforming gives 
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and Fourier 
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where is the area of the surface unit cell, and {G} is the set of reciprocal lattice vectors of the overlayer structure. 

Equation ( p^ ) demonstrates that a coherent excitation of an adsorbate overlayer cannot couple to all resonant 
substrate modes, in contrast to the case of an isolated adsorbate.B Rather, the lattice imposes a selection rule. 
Interference among the adsorbates in the overlayer prevents energy from being radiated into any substrate mode 
whose in-plane wavevector component is not commensurate with the adsorbate reciprocal lattice. Furthermore, the 
factor of 1/Ac in Eq. (|l^) shows that dense adlayers drive the substrate more strongly than diffuse adlayers, or, in 
the extreme limit, than an isolated adsorbate. Thus the coupling is more efficient for a dense overlayer and leads to 
a faster damping, as shown below. 

Combining Eqs. (|^) and (|lo|), inverse Fourier transforming, and evaluating the result at x = -R^ gives a simple 
expression directly relating the substrate and adsorbate degrees of freedom: 



Ux{Ri,ljj) = Si{uo) - Ux{Ri.,Uj) 
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Dxx{q + G,uj) 



(11) 



The response function, x('fj^)j defined as the reciprocal-space version of the function appearing in Eq. (||), can 
be expressed in terms of the substrate Green's function, Dxx{q,i^), and fundamental parameters of the system by 
combining Eqs. (0) and (O), and rearranging the result to be in the form of Eq. (0). Thus, 



where 



1 + T{q,cu) 



l + Riq,Lu) 



{l + R{q,u;)y+P{q,oj) 



lUJQ 



mujal{q,u;) 



{l + Riq,oj)y+P{q,cj) 
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and R{q,uj) and I{q, uj) are the real and imaginary parts, respectively, of T{q,uj). The isolatcd-adsorbate extreme 
is obtained from Eq. (O) in the limit of infinitely long overlayer lattice vectors. The Brillouin zone area, (27r)^/Ac, 
becomes an infinitesimal and the sum over reciprocal lattice vectors becomes an integrahc 



raujr. 
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The imaginary part of x((f, u), Im[x(<f, i-i^)], is the measurable absorption spectrum. If we assume that T{q,Lo) is 
sufficiently smooth and that the damping rate, 7, is much smaller than the resonant frequency, then Im[x(9, i^)] is 
a sharply peaked function centered at the renormalized frequency (Do (the experimentally measured frequency) and 
having characteristic width 7. Under these conditions, we can approximate uj ^ wq, which gives x{q^^) the form of a 
damped harmonic oscillator: 



x((f, cj) « \—muj^ + mujQ — iujmj^ ^ , 



(15) 
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where 



l + R{q,UJo) 



(1 
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R{q,LLio)) - 
uJoIiq, ujq) 



{l + R{q,u;o)r+P{q,u;o) 



(16) 
(17) 



For the remainder of this paper we assume that dispersion in the plane is neghgible. This assumption is governed by 
the length scale for inter-adsorbate interactions, and is expected to be valid for all but the highest surface coverages. 
Thus all analyses are made for the q = limit. This restriction could be lifted without difficulty, if warranted. 

Determining the dynamical response of the adlayer has been reduced to solving for the Green's function of a flat, 
semi-infinite, isotropic elastic medium. This has been carried out in great detail by Maradudin, et al.u Evaluating 



their expressions for D^^ at z 
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[{a - (3)k^ + (a + f3)p- 



(19) 
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and a — p{cj^ — ) and f3 ~ pc^. Equations (|l8|)-(p0|) use the same naming conventions as in Appendix B of PR.0 The 
analysis that has been presented so far for x-polarized adsorbate motion also holds for z-polarized motion, provided 
that Uz and D^z replace Ux and D^x, respectively, throughout, where Dzz is given by 



Dzzik,Lu) = -i 



-Pt 



2g(fc2,w) 



(21) 



According to Eq. (|T^), Dxx{k,uj) is to be evaluated only at reciprocal lattice vectors, k = G. This function is 
pure imaginary unless pt or pj^ is imaginary. Consider the shortest non-zero reciprocal lattice vector, Gq. Assuming 
a square adsorbate lattice, jGol = 27r/a, where a is the adsorbate lattice constant. The condition for pt.l to be 
imaginary is lo/ct^l < At Gq, this is equivalent to 



a < Xt,l , 



(22) 



where \t,l = '^t^ct.l/^ is the wavelength of a substrate phonon at or near the resonant frequency. However, since 
the elastic-continuum description of the substrate is only valid at the long-wavelength end of the phonon spectrum, 
Eq. (^2|) will be true for the present theory unless the adsorbate lattice is extremely dilute. Suppose, for example, 
that A = lOcbuik: where Obuik is the lattice constant of the bulk substrate material. Since surface coverage is given 
hy — abuik/*^^' l^'l- ( p2| ) will be violated only for coverages less than 1%, which is below the resolution limit of 
current surface-sensitive measurement techniques. Thus, for all systems presently accessible to experiment, all non- 
zero reciprocal lattice vectors contribute only to the real part of T(uj), and only G — contributes to the imaginary 
(i.e., damping) part. We call the regime specified by Eq. ( p^ ) the "dense adlayer" limit, and it is given in terms of 
surface coverage by 
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where lod = ct,l 7r/abuik is the substrate Debye frequency. 
In the dense adlayer limit, the imaginary part of r((jJo) is 

/(c.o) = ^i?..(0,c.o) 



mujo 

AcpCT 



(24) 
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TABLE I. Empirical physical parameters for a half-monolayer of carbon monoxide on the copper (100) surface. 



Quantity 



Value 



^0 (10^^ s-i) 
m (lO^^^* g) 
a (10~* cm) 
A, = (10-1' 

P (g/cm^) 
CT (W^ cm/s) 
CL (10^ cm/s) 



6.03 
4.65 
3.61 ^ 
1.30 
8.93 ^ 
2.91 
4.34 " 



'■Reference 7 



"Reference 



'^Derived from elastic constants in Reference I 



which, when inserted in Eq. (p7|), gives a damping rate of 

mUJn 1 
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AcpCT {1 + R{ujo)y + P{uJo) AcpCT 



1 + i?(c^o) 



AcpCT 



(25) 



where the expression for lDq in Eq. ( |16D has been used. Typically, |i?(aJo)| << Ij so that the bracketed term in Eq. ( |25| ) 
can be neglected. Thus the vibrational damping rate of a dense adlayer can be simply computed from empirically 
available information about the system. Note that the expression for 7 in Eq. (25) is strongly dependent on the 
surface coverage through the parameters Ac, wq, and R{ujo). 

The real part of T(u;o) is a sum over all non-zero reciprocal lattice vectors. For real systems, the substrate is not a 
continuum, but rather a discrete lattice. This places an upper bound on the summation limit for the computation of 
R{ujo): the sum includes all reciprocal lattice vectors of the overlayer structure that lie within the first BriJlouin zone 
of the bare surface. Surface wavevectors outside this region do not correspond to new substrate phonons.EI 

The Green's function DxxiG,ujQ), for G ^ 0, can be simplified in the dense-adlayer limit by expanding to lowest 
order in {ujq / cl ,t G)"^ , where G = \G\: 



Dxx{G,ujo) 



1 



pc^G 



Gl Gl 
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G2 G2 2(1-(ct/cl)2) 



Since the term in brackets is of order unity, R{ujo) is approximately given by 



mujQ 
Acpc^G 



(26) 



(27) 



ff the overlayer unit cell is small (i.e., high coverage), then the sum will contain only a few reciprocal lattice vectors 
before reaching the upper bound set by the substrate lattice. Thus, since ojo/ctG ^ 1, i?(wo) will be small. As the 
unit cell gets larger, more terms are included in the sum. However, the factor 1/Ac decreases, keeping i?(ijJo) small. 

To test the validity of this theory, we now apply it to a system for which there are experimental data: carbon 
monoxide on the (100) surface of copper at half-monolayer coverage. The in-plane frustrated translational mode has a 
measured frequency of 32 cm"^ , well within the range of long- wavelength copper phonons. Recently, the lifetime of this 
mode was measured to be 2.3±0.4 ps,Ll which corresponds to a damping rate (0.43±0.07) x 10^^ s"^. Other parameters 
for this system are compiled in Tabic |. Using these data, Eq. (p5|), with the bracketed term neglected, predicts the 
frustrated translational damping rate to be 7 = 0.50 x 10^^ s~^, in excellent agreement with the experimental value. 

The bracketed term in Eq^23) produces a small downward correction. To show this we now compute R{uJo) 
using the expression in Eq. (^). The CO overlayer forms a square lattice with c(2 x 2) symmetry relative to the 
underlying copper (100) surface. Therefore, the only non-zero adsorbate reciprocal lattice vectors that lie within the 
first Brillouin zone of the bare surface are [±27r/a, 0] and [0, ±27r/a], leading to a sum with only four terms. The 
resulting expression for R{ujo) for this system is 



2Trpc^a \ ' 1 — [ct/clY' 



= 0.0376 



(28) 



Using this value, the corrected damping rate is 0.48 x 10^^ s ^ . Furthermore, the estimated frequency renormalization 
factor is 0.956, which, as anticipated, is close to unity. Note that the value of R{uJo) was computed using the measured 
frequency and not the bare frequency. However, since the renormahzation is small, this approximation is acceptable. 
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In prevkjus work, we performed an atomistic investigation of this system treating the interactions quantum me- 
chanicallyE The computed frustrated translational frequency and damping rate were 27 cm~^ and 0.33 x 10^^ s~^, 
respectively. Thus the present elastic continuum theory, the atomistic theory, and the experiments all give a consis- 
tent, mutually reinforcing picture of the low-frequency adsorbate vibrational relaxation dynamics; namely, a layer of 
oscillating adsorbates coupling collectively to substrate phonons at the resonant frequency. 

While much research on adsorbed surfaces has focused on ordered overlayer structures, it is also important to 
consider how the generic case of a disordered or defective overlayer affects the damping rate formula of Eq. (^5|) . For 
a dense adlayer (as defined above), the relaxation is governed by atomic motions in which all atoms of a given layer 
oscillate in phase, for both the adlayer and the substrate. Thus the relevant behavior is one-dimensional along the 
surface-normal direction, and the details of the in-plane structure are not resolved. As a result, the damping rate 
formula, including its coverage and frequency dependence, is unaffected by disorder in the dense-adlayer limit. In this 
analysis we have assumed that the adsorbates are uniformly distributed on the surface (i.e., no islands). If islands 
form, then the local adsorbate density is higher than the average density, and we expect the resonant damping rate 
to reflect the former value. A more thorough, quantitative treatment of disorder in this problem will be presented in 
a forthcoming article. 

Resonant coupling between an adsorbate overlayer and an elastic substrate provides a damping force for the decay 
of low-frequency vibrations. We have presented a classical, analytical theory of this mechanism which predicts 
vibrational damping rates in excellent quantitative agreement with measured values for systems that have been 
studied experimentally. 

Furthermore, we have shown that, even for the smallest measurable coverages, vibrating adsorbates behave collec- 
tively when coupling to the substrate. Due to interference effects, a layer of adsorbates oscillating in phase couples 
only to resonant phonon modes propagating perpendicular to the surface. Only when neighboring adsorbates are 
separated by distances larger than the wavelength of resonant phonons will new propagation directions be included. 
The isolated-adsorbate limit, treated by Persson and Ryberg,0 becomes valid only when many resonant-phonon wave- 
lengths fit between neighboring adsorbate sites. This limit is well below the measurable threshold. Finally, we have 
also shown that the coupling is strongly coverage dependent. A dense overlayer drives the substrate more efficiently 
than does a dilute overlayer. 

We thank J. P. Toennies for useful discussions during the early stages of this work. Financial support for this 
project was provided by the National Science Foundation under grants DMR 97-02514 and DMR 93-13047, and by 
the Petroleum Research Fund of the American Chemical Society under grant 32007-G5. 
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